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1 Introduction 



As is well known, one-dimensional (ID) quantum integrable spin chains can be broadly 
classified into two types depending on their range of interaction. One such class consists 
of spin chains having only local interactions like nearest or next-to-nearest neighbour 



model etc. are examples of this type of spin models with short range interaction [1-3]. 
However, there exist another class of ID quantum integrable spin systems for which 
all constituent spins mutually interact with each other through long range forces [1, 
4]. Haldane-Shastry (HS) spin chain [5,6], Polychronakos spin chain (also known as 
Polychronakos-Frahm (PF) spin chain) [7-9] and Frahm-Inozemstsev (FI) spin chain 
chain [10, 11] are well known examples of quantum integrable spin systems with long 
range interaction. This type of spin chains have attracted a lot of attention in recent 
years due to their exact solvability and applicability in a wide range of subjects like 
generalized exclusion statistics and fractional quantum Hall effect [12-15], SUSY Yang- 
Mills theory and string theory [16-20], Yangian quantum group [21-28] etc. For the 
simplest case of root system, the Hamiltonians of su(m) invariant, ferromagnetic 

type HS, PF and FI spin chains can be written in a general form like [4-11] 



Here denotes the position of the j-th lattice site, spins can take m possible values on 
each lattice site, Pjk represents the exchange operator which interchanges the spins on 
the j-th and fc-th lattice sites and the strength of this exchange interaction is determined 
by the function /(£,■ — It should be noted that, while the lattice sites of HS spin 
chain are uniformly distributed on a circle, the lattice sites of PF and FI spin chain 
are non-uniformly distributed on a straight line. Supersymmetric generalizations of the 
spin chains (1.1) have also been studied in the literature [13,29-31]. 

The difference between quantum integrable models with short range interaction and 
long range interaction manifests very prominently through their connection with classical 
vertex models. It is well known that the Hamiltonians of ID quantum integrable spin 
chains with short range interaction are related to the transfer matrices of some exactly 
solvable 2D vertex models in statistical mechanics [32,33]. However, there exists no 
simple relation between the energy levels of such ID spin chains and the energy functions 
of corresponding 2D vertex models. On the other hand, a direct connection between 
the energy levels of ID quantum integrable spin chains with long range interaction 
and energy functions of ID vertex models can be inferred from some works in the 
literature [24,26,34]. Indeed, by applying a rather general technique based on Yangian 
quantum group symmetry, it has been shown recently that the energy levels for a class 
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of HS like spin chains as well as their supersymmetric generalizations exactly coincide 
with the energy functions of some ID inhomogeneous vertex models in classical statistical 
mechanics [35]. 

The reason behind the close connection between ID quantum spin chains with long 
range interaction and ID classical vertex models can be traced from the underlying 
Yangian quantum group symmetry of the spin Hamiltonians, which allows one to write 
down their exact spectra in closed forms even for finite number of lattice sites. In 
the following, we shall briefly recapitulate the consequences of Yangian symmetry for 
the special case of nonsupersymmetric su(m) ferromagnetic spin chains. A class of 
irreducible representations of the Y(sl m ) Yangian algebra, known as 'border strips' or 
'motifs', completely span the Hilbert space of some Yangian invariant spin systems with 
long range interaction [4, 13,21-28]. For the case of a spin system with N number of 
lattice sites, boarder strips are denoted by (ki, ki, ■ ■ ■ , k r ), where fcj's can be taken as any 
positive integers satisfying the relation YH=i ^* = N. Thus the Hilbert space associated 
with a Yangian invariant spin system may be decomposed as 

V= J2 ® V(kuk 2 ,..,k r) , (1-2) 
k 1 +k 2 +...+k r =N 

where V^ lt k 2 ,...,k r ) denotes the irreducible vector space associated with boarder strip 
(&4, &2, . . . , k r ). Let us restrict our attention to a class of Yangian invariant spin systems, 
for which the complete set of energy levels can be written in the following way: 

r-1 

E( ku k 2 ,...,k r ) = y^^JvQ')|j=jfi , (1-3) 
1=1 

where Ki = Y^i=i and £n{j) is an arbitrary function of the discrete variables j and N. 
Due to Yangian symmetry of this class of spin systems, the multiplicity of the eigenvalue 
E(k 1 ,k 2 ,-,k r ) i n Eq. (1.3) coincides with the dimensionality of the vector space V(k u k 2 ,...,k r )- 
Hence, this type of energy eigenvalues are highly degenerate in general. It should be 
noted that, all energy levels of su(m) ferromagnetic HS, PF and FI spin chains can 
be generated through Eq. (1.3) as special cases, where 5jv(j)' s are given by 'dispersion 
relations' of the form [21-28,34] 

( j(N-j), for the HS chain 

E N (j) = 1 j, for the PF chain (1.4) 

[ j( c + j ~ 1)> w ^h c > 0, for the FI chain. 

The spectra of all Yangian invariant spin chains associated with Eq. (1.3) can be 
written in an elegant alternative form by using the energy functions of corresponding ID 
inhomogeneous vertex models [35]. Let us consider a class of ID classical vertex models 
with (JV+1) number of vertices, which are connected through iV number of intermediate 
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bonds. Each of these bonds can take one of the m possible states. Therefore, any state 
for such vertex models can be represented by a path configuration of finite length: 

s = {s 1 ,s 2 ,...,s N }, (1.5) 

where G {1, 2, • • • , m} denotes the spin state of the i-th bond. For constructing vertex 
models corresponding to Yangian invariant spin chains whose energy levels are given by 
Eq. (1.3), let us define an 'energy function' related to the spin path configuration s as 

N-l 

E(s) = J2^n(j)S( Sj ,s j+1 ), (1.6) 

where 5(sj, Sj+i) is taken as 

1*Z> I (L7) 

It should be noted that, the same function £n(j) has appeared in Eqs. (1.6) and (1.3). 
It can be shown that, for all possible functional forms of such £nU), the spectrum of the 
su(m) ferromagnetic spin chain generated by eigenvalues (1.3) (along with the degen- 
eracy factors of these energy levels) completely matches with that of the corresponding 
vertex model generated by energy functions (1.6) [35]. In particular, the full spectra 
of su{m) ferromagnetic HS, PF and FI spin chains can be reproduced from the energy 
function (1.6), by substituting to it the corresponding forms of dispersion relations given 
in Eq. (1.4). Hence, it is possible to extract important informations about the statisti- 
cal properties of the spectra of these spin chains by analyzing the distribution of energy 
functions for the related ID vertex models. With the help of numerical studies, it has 
been conjectured for a long time that the energy level densities for this type of spin 
chains asymptotically follow the Gaussian distribution as the number of lattice sites 
become very large [11,36-41]. By using the above mentioned connection with ID vertex 
models, very recently this conjecture has been proved analytically for the case of su(m) 
HS, PF and FI spin chains associated with the A N _i root system [34]. 

In this context, it is interesting to ask whether the Gaussian level density distribution 
of HS, PF and FI spin chains is simply a consequence of their Yangian quantum group 
symmetry or this result crucially depends on the specific forms of the dispersion relations 
given in Eq. (1.4). To answer this question, in this article our aim is to explore the level 
density distribution of the energy function (1.6) for a large class of ID vertex models, 
for which En{j) can be chosen as an arbitrary polynomial of the variables j and N. 
Even though it may not be possible in general to express the Hamiltonians of the spin 
chains corresponding to such vertex models in a simple form like (1.1), one can formally 
express those Hamiltonians in dyadic notation through the basis vectors of the Hilbert 
space V in the following way. Let us denote the basis vectors of the subspace V(ki,k2,...,kr) 
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in Eq. (1.2) with the notation |V , fci,fc 2 ,-,fcr;a)> where a G {1,2,--- ,pk}, and pk is the 
dimension of V(fci,fc 2 ,...,fcr)- The Hamiltonians for the above mentioned class of Yangian 
invariant spin chains may now be written as 

W = ^2 ff(fci,fc 2 ,...,fcr.) ^ !V'fci,fc2,...,fer;a)(V'fci,fe2,...,fc r ;al > (I- 8 ) 

fei+fe2+...+fcr=A f a=l 

where E^^,...,/^) i s given by Eq. (1.3), with En{j) being an arbitrary polynomial func- 
tion of j and N. Hence, in the following, our aim is also to explore the level density 
distribution for the class of spin chains represented by % in Eq. (1.8). 

The arrangement of this article is as follows. In Sec. 2 we briefly recapitulate the 
transfer matrix approach of Ref . [34] , where it has been found that the level densities for 
the case of nonsupersymmetric HS, PF and FI spin chains and related ID vertex models 
asymptotically follow the Gaussian distribution as the size of these systems become very 
large. In Sec. 3 we show that the above mentioned approach can be applied to analytically 
compute the level density distribution for a wide range of ID vertex models and related 
spin chains, for which the dispersion relations are taken as arbitrary polynomials of the 
variables j and N. In this section, we also find out the order of mean and variance 
for the energy functions or energy levels associated with such dispersion relations. In 
Sec. 4 we numerically study the level density distribution for some ID vertex models with 
different choice of the related parameters and find that these numerical results are in 
complete agreement with the analytical prediction of the previous section. In Sec. 5 we 
summarize our results and make some concluding remarks on the possible connection 
between Yangian quantum group symmetry of a spin chain and the asymptotic form of 
its level density distribution. 



2 Level density for known HS like spin chains 

In this section, we review the transfer matrix approach for calculating the level density 
distributions in the cases of ferromagnetic su(m) HS, PF and FI spin chains related to 
the An -i root system [34]. A key role in this approach is played by the normalised 
characteristic function of the level density distribution. Utilizing the equivalence of 
these spin chains with ID vertex models, such normalised characteristic function may 
be defined as 

«t) = (e^>, (2.1) 

where E(s) is the energy function (1.6) of the vertex model associated with the spin 
chain, p and a denote the mean and the variance of this energy function respectively, 
and the notation (O) represents the average of O over all spin path configurations. If 
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this characteristic function satisfies the asymptotic relation 



lim 0jv(t) = e" 

TV-i-oo 



(2.2) 



then it can be shown that [42] at iV — > oo limit the normalised level density approaches 
to a Gaussian distribution of the form 



G(E) = ^=e 



Therefore, instead of directly working with the level density of a spin chain, one can 
study the N — > oo limit of the related characteristic function. Inserting the form of 
E(s) given in Eq. (1.6) to Eq. (2.1), it is easy to express 0jv(t) as 

} N (t)=m- N e-^Z N (e'), (2.4) 

where the 'partition function' Zjv(e^) is defined through a product of transfer matrices 



as 



Z N ^)=m N -' J2 [TMT(u 2 )...T(u N ^)] nnl 

n,n'=l 

Here the transfer matrix T(ujj) is given by T(tUj) = T(oj)\ u=Uj , with 



(2.5) 



T{u) 



m 



ojj = e tljt and 



7i 



(1 u m ■ 

\l ... . 

SnU) 



(2.6) 



1 / 



ma 



(2.7) 



Since this transfer matrix can be diagonalised through an unitary transformation, it 
may be written as 

T( Uj ) = U( Uj )D( Uj )lPM, (2.8) 
where U(ujj) is an unitary matrix with elements given by 



and D(u)j) is a diagonal matrix: 

D(wj) = diag (Ai(wj-), • • • , A m (u^)) , 

with 



(2.9) 



(2.10) 



(2.11) 
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With the help of Eqs. (2.4), (2.5) and (2.8), one can express 07v(t) as 

-%iit i a m 

h{t) = Yl M nn>(t), (2.12) 

n,n =1 

where the matrix M is given by 

M(t) = U M D M B 1 (t)---D ( Wjv _ 2 ) B N _ 2 {t)D (ou^) (c^v-i) , (2.13) 

with 

B j (t) = rf(u, j )U(«, j+1 ). (2.14) 

It should be observed that, Eq. (2.12) is derived without assuming any specific form of 
the dispersion relation. Furthermore, for any functional form of Sjy(j), the mean and 
the variance of the energy function E(s) in Eq. (1.6) can be expressed as [34] 



1 / ^\ N - 1 



and 

, iV-l N-2 

a* ee ((E(s) - rf) =U-—\-J^ E N {jf - - SnWnU + 1) 



(2.16) 



However, for determining the N — > oo limit of 0jv(£) in Eq. (2.12), it is necessary 
to choose some specific form of the dispersion relation. Inserting the specific forms of 
dispersion relations given in Eq. (1.4) to Eqs. (2.15) and (2.16), one can explicitly obtain 
pi and cr 2 (as some functions of m and N) for the cases of HS, PF and FI spin chains. 
By plugging such explicit forms of a in Eq. (2.7), it can be shown that the asymptotic 
behaviours of jj and jj — 7 J+ i for large values of N are given by the following relations 
in the cases of all these three spin chains [34]: 

15 = 0{N-W), (2.17) 

and 

U ~ 7;+i = 0(N~^) , (2.18) 
where j e {1, 2, . . . , TV - 1} and the notation f(N) = 0(N a ) implies that 

lim \N~ a f(N)\ =t, (2.19) 

with r being some non-negative real number. With the help of Eqs. (2.9) and (2.17), it 
is easy to show that 

U(uj) = R + 0{N^ 2 ) , (2.20) 
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where R is a constant unitary matrix with elements given by 

R nn , = _L e -2Tinn'/m _ ^.21) 



By using Eqs. (2.9), (2.14) and (2.18), one can also determine the asymptotic form of 
Bj(t) as 

Bj(t) = l + 0(N- 3/2 ). (2.22) 

Substituting the relations (2.20) and (2.22) to Eq. (2.13) and taking the N — > oo limit, 
one obtains 

lim M(t) = R\(t)R\ (2.23) 

JV->oo 

where A(t) is a diagonal matrix given by 

A(t) = diag(A 1 (t),...,A m (t)), (2.24) 

with 

JV-l 

A fc (t) = lim FT \ k (u>j) . (2.25) 

JV— >-oo 

Taking the N ->■ oo limit of Eq. (2.12), and subsequently using Eqs. (2.23), (2.21) and 
(2.25), it is easy to find that 

JV-l 

lim <Mt) = e-^A m {t) = e~¥ lim TT \ m {uj) . (2.26) 

JV— >oo JV->oo - LJ - 

By putting k = m in Eq. (2.11), one obtains X m (ujj) = — YlT=o e * r<7j ' f rom which it 
is clear that \ m {uij) — > 1 in the limit jj — > 0. Expanding In \ m (u>j) in a power series 
of 7^, and subsequently using the expressions (2.7), (2.15) and (2.16) along with the 
asymptotic forms (2.17) and (2.18), one gets the relation 



'JV-l 



\n{l[\ m (^)\ = ^-^+0(N-i), (2.27) 

3=1 



which in turn leads to 

JV-l 

ifit t 



lim TT \ m (wj) = e^~^. (2.28) 



JV^OO 

J'=l 



Substituting this result to Eq. (2.26), one finally obtains the desired asymptotic form 
(2.2) for the normalised characteristic function. Thus it is established that the level 
density distribution for su(m) ferromagnetic HS, PF and FI chains asymptotically follow 
the Gaussian pattern (2.3) for large number of lattice sites. 
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3 Level density for vertex models with polynomial 
type dispersion relations 



From the discussions of the previous section it is evident that the asymptotic forms 
of 7j and jj — 7j+i, given by Eqs. (2.17) and (2.18) respectively, play an important 
role in determining the level density distribution for the cases of HS, PF and FI spin 
chains. In this context it is interesting to observe that, even though jj in Eq. (2.7) 
explicitly depends on the functional form of £/v(i), the asymptotic form of jj (or, jj — 
7j+i) is completely same for all three dispersion relations appearing in Eq. (1.4). This 
observation suggests that some universality might be present in these asymptotic forms 
of 7j and 7j — 7j+i- As a result, the transfer matrix approach might be applicable for 
finding out the level density distribution of a much larger class of vertex models and 
associated quantum spin chains. 

Indeed, in the following our aim is to consider a large class of vertex models where 
the dispersion relations can be taken as arbitrary polynomial functions of the variables 
j and N, and show that the asymptotic forms of jj and jj — 7^+1 are universally given 
by Eqs. (2.17) and (2.18) respectively for all such cases. To begin with, we write down 
the general form of such polynomial type dispersion relation as 

r s 

Mj) = ££/«,* (3-i) 

a=0 /3=0 

where r and s are some non-negative integers, and f a ,^s are arbitrary real numbers. 
For the sake of uniquely expressing a polynomial type dispersion relation in the form 
(3.1) with minimum possible values of r and s, we assume that (without any loss of 
generality) there exist at least one value of /3 e {0, . . . , s} such that f r $ 7^ and at 
least one value of a G {0, . . . , r} such that f ayS 7^ 0. The degree of the two-variable 
polynomial £n{j) in Eq. (3.1) is denoted by x, which takes an integer value within the 
range: ^ x ^ r + s - Evidently, /a^'s appearing in Eq. (3.1) satisfy the conditions 

faj = 0, for a + /3 > x, 

7^ 0, for at least one pair of a, (5 with a + (5 = x ■ (3-2) 

By using these conditions, it is easy to show that 

X ^ max{r, s}. (3.3) 

It may be noted that the general form of dispersion relation (3.1) yields all dispersion 
relations given in Eq. (1.4) as some special cases. For example, the choice r = 0, s = 1 
along with / ,o — 0, /o,i — 1 gives the dispersion relation for PF spin chain and the 
choice r = 1, s = 2 along with / ,o = /o,i = /i,o = /i,2 = 0, / 0>2 = ~fi,i = -1 gives the 
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dispersion relation for HS spin chain. Clearly, x = 1 for the vertex model corresponding 
to PF spin chain and x — 2 for the vertex models corresponding to HS and FI spin 
chains. 

Even though it is easy to find out the explicit form of a 2 by using Eq. (2.16) for 
the above mentioned type of special cases associated with small values of r and s, it is 
quite cumbersome to do so for the dispersion relation (3.1) with high values of r and s. 
Therefore, instead of trying to compute the exact form of a 2 explicitly, at present our 
strategy is to estimate the order of a 2 for the general form of dispersion relation (3.1). 
To this end, we consider the 'big notation' for the order of a function: f(N) = Q(N a ), 
which implies that 

Jim \N- a f(N)\=r, (3.4) 

where r is a positive real number. It should be observed that, in contrast to the case of 
the 'big O' notation defined through Eq. (2.19), at present r is not allowed to take the 
zero value. As a result, the expression f(N) = Q(N a ) carries more precise information 
about the asymptotic form of f(N) in comparison with the expression f(N) = 0(N a ). 
Indeed, while f(N) = 0(N a ) implies that f(N) is bounded above by iV a (up to a 
constant factor) asymptotically, f(N) = Q(N a ) implies that f(N) is bounded above 
and below by N a (up to constant factors) asymptotically [43-45]. In the following, we 
shall use both of these notations for the order of a function according to our convenience. 

At first, let us examine the order of £n{j) in two possible cases: i) j = 0(1), where 
iV — > oo limit is taken by keeping a fixed value of j, and ii) j = Q(N), where N — > oo 
limit is taken in such a way that j/N tends to a fixed value. Since there exists at least 
one value of j3 e {0, . . . , s} for which f r $ ^ in (3.1), we easily obtain the order of 
£n{j) for case i) as 

£ N {j) = Q(N r ) . (3.5) 
Next we consider the case ii), which is defined through the limiting procedure 

Jim T7 = P, (3-6) 

where p is a real parameter taking value within the range: < p ^ 1. Due to Eq. (3.6), 
one may write £n{j) ~ £n(pN) for large values of N. Substituting pN in the place of 
j in Eq. (3.1), we obtain 

r s 

s N (pN) = j2Y,f^p" Na+ "- ( 3 - 7 ) 

a=0 /3=0 

Let us define a new variable as £ = a + /3. Making a change of the summation variables 
in Eq. (3.7), it can be expressed as 

r+s min{£,s} 

£n( P n) = J2 E h-fut**- ( 3 - 8 ) 

(=0 /3=max{0£-r} 
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Collecting the coefficients of in the above equation, we get 



r+s 



£N(j>N) = Y,9e(p)N*, (3-9) 
where g${p) is a polynomial of p given by 

min{t;,s} 

9dp)= E A-A/J^- (3-10) 

/3=max{0,5— r} 

Since from Eq. (3.2) it follows that g^(p) = for £ > x, one can rewrite Eq. (3.9) as 

x 

S N ( P N) = Y,9 6 (P)^. (3.11) 

£=o 

With the help of Eqs. (3.6) and (3.11), we obtain 

lim \N-*£ N (j)\ = lim \N~* £ N (pN)\ = \g x (p)\ . (3.12) 

iv— »oo N— >oo 

Since g x (p) is a polynomial of p, it is bounded within the range < p ^ 1. If g x (p) 7^ 
for some value of p, then by using Eqs. (3.4) and (3.12) we obtain the order of £n{j) as 

£ N (j) = 9(iV*). (3.13) 

On the other hand, if g x (p) = for some value of p, then Eq. (3.11) reduces to £ N (pN) = 
Y%=o9t{p)N^, wnere X < X an d fi'x(p) 7^ 0- Hence, for the case g x (p) = 0, we obtain 
the order of £n{j) as 

£ N (j) = e(N*), (3-14) 

where x < X- Replacing £ by x in Eq- (3.10) and using Eq. (3.3), one can get an explicit 
expression for g x (p) as 

s 

9xip) = E fx-MfP- (3-15) 

/3=x-r 

Since <? x (p) is a polynomial in p of degree less than or equal to s, it can have at most s 
number of distinct zeros within the interval < p ^ 1. For such zero points of g x (p), 
the order of £n{j) is evidently given by Eq. (3.14). However, the zero points of g x (p) 
clearly form a set of zero measure within the range < p ^ 1, and, therefore, the order 
of £nU) is given by Eq. (3.13) for generic values of p. 

Let us now compare the orders of £n{]) given in Eqs. (3.5), (3.13) and (3.14) for 
various possible cases. Since due to Eq. (3.3) it follows that x > r ? the order of £n{j) 
given by (3.13) is clearly the highest one among all of these cases. Therefore, at N — > 00 
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limit, the expression of a 2 in Eq. (2.16) would get dominant contribution from those 
values of j which satisfy the conditions j = @(N), g x (p) ^ and lead to Eq. (3.13). 

Next we explore the stability condition for the order of £n{j) given in Eq. (3.13), 
under the variation of parameter j (or, equivalently, p). Since the polynomial g x (p) can 
have at most s number of distinct zeros within the interval < p ^ 1, we can always find 
two distinct numbers p\ and p 2 within this interval such that g x (p) ^ for p 1 ^ p ^ p 2 
and 

p2-pi = —r • (3.16) 
s + 1 

Due to Eq. (3.6), one can write N pi j\ and Np 2 ~ j 2 for large values of N, where 
ji, j 2 are some integers. Since g x (p) ^ for p 1 ^ p ^ p 2 , it is clear that Eq. (3.13) is 
satisfied for any j e {ji, ji + 1, . . . , j 2 }. Thus, we can write 

S N {ji+Aj)=Q{N^), (3.17) 

where Aj may take any value within the range j ^ j 2 — j\. By using (3.16), we 
obtain j 2 — j 1 N(p 2 — p±) — which yields the order of j 2 — ji as 

j2-ji = Q(N). (3.18) 

Eqs. (3.17) and (3.18) imply that, there exists at least one value of j (denoted by j\ 
here) for which £n{j) = Q(N X ) and this asymptotic behaviour of £n{j) is preserved for 
a variation of j of order N . This type of stability in the asymptotic behaviour of SnU) 
will play an important role in our analysis in finding out the order of a 2 . 

For the purpose of finding out the asymptotic behaviour of a 2 , it is also needed to 
calculate the order of £n{j + 1) — £n{j)- By using Eq. (3.1), we obtain 

r s /3-1 

£' N (j) = £ N (j + 1) - £ N (j) = Pc »f°>e Na ? ' 

a=0 13=1 v=0 

where ^C u = vl ^!_ v y m ■ Interchanging the order of summations over (3 and v in the above 
equation, and renaming these summation variables, we express £' n{j) in exactly same 
form as £n(j) in Eq. (3.1): 

a=0 (3=0 

where r' < r, s' = s — 1 and f' a p is given by 

s 

f'^= E f^'Cp. (3.20) 

v=P+l 
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Consequently, our analysis for finding out the order of £n{j) w °uld also be applicable 
for the case of £' N (j). Let us denote the degree of £' n{j) in Eq. (3.19) by x' ■ With the 
help of Eq. (3.20) and the analog of Eq. (3.2) for the present case, it is easy to find that 

X'<X-1- (3.21) 

In analogy with the polynomial g x (p) in Eq. (3.15), one can define the polynomial g x >(p) 
for the present case as 

s' 

9AP)= E fr-wPP- ( 3 - 22 ) 

0= x '-r> 

Proceeding in the same way as we have done to find out the maximum possible order of 
£n{j) i n Eq. (3.13), it is easy to show that the order of £' n{]) becomes maximum when 
j satisfies the relations j = Q(N) and g x >(p) ^ 0, and this order of £' n{j) is given by 

£' N (j)=e(Nx'). (3.23) 

Next, to calculate the order of a 2 for any dispersion relation of the form (3.1), we 
rewrite Eq. (2.16) as 

a 2 = h + I 2 + J 3 , (3.24) 



where 



i / i \ N - 1 
'i=12 i-^lEMi) 2 , 0.25a) 



i=i 

N-2 



If 1 \ 

h = -gf 1 - —A J2^U)£'nU), (3.25b) 

3=1 

h = Ul-^)£N{N-l) 2 . (3.25c) 

Since terms like £n{j) 2 appearing in the r.h.s of (3.25a) are positive for all values of 
j, their contribution in Ii can not cancel out each other. Therefore, by applying Eqs. 
(3.17) and (3.18), we find that that h = 9(7V 2 * +1 ). Due to Eqs. (3.21) and (3.23), it 
is evident that the order of I2 i n (3.25b) must be less than that of I±. Moreover, since 
the expression of 1% in (3.25c) does not contain any summation variable, the order of it 
must be less than that of I\. Consequently, the dominant contribution in the r.h.s. of 
Eq. (3.24) comes from the term I\. Thus, we obtain the order of a 2 as 

a 2 = Q(N 2x+1 ). (3.26) 

Since the parameter r in Eq. (3.4) takes a nonzero value, it is possible to 'invert' the 
relation (3.26) and calculate the order of 1/er as 

I = 0(jV-*-5). (3.27) 

<7 
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Now, we are in a position to find out the orders of 7,- and 7j — 7j+i for any polynomial 
type dispersion relation of the form (3.1). Let us first consider the case for which 
j satisfies the relations: j = Q(N), g x (p) 7^ 0, and the order of £n{j) is given by 
Eq. (3.13). By using Eqs. (2.7), (3.13) and (3.27), we obtain the order of 7, for this case 
as 

7i = e(j\r*). (3.28) 

Since the order of £n{j) becomes maximum for the above mentioned case, it is evident 
that the order of 7^ can not exceed the value given in the r.h.s of Eq. (3.28) for any 
possible choice of j. In other words, the asymptotic form of 7, at N — > 00 limit can not 
exceed (up to some multiplicative factor) for any possible choice of j. Hence, by 

using the 'big O' ordering defined in Eq. (2.19), we obtain a relation valid for any value 
of j as 

7j = O(TV-I) . (3.29) 

Next, for calculating the order of 7j — 7^+1, we proceed in exactly same way and consider 
the case for which j satisfies the relations: j = Q(N), g x '(p) 7^ 0, and leads to the 
maximum order of £'n{j) given by Eq. (3.23). Using Eqs. (2.7), (3.23) and (3.27), we 
obtain the order of 7,- — t,- +1 for this case as 

7i - 7i+1 = 6(^-^-4). (3.30) 

In analogy with Eq. (3.28), the above equation implies that the asymptotic form of 7, — 
7j + i can not exceed jV x '~ x ~2 for any possible choice of j. Moreover, due to Eq. (3.21), 
it follows that x' — x — | < — §. As a result, we can express the order of 7, — 7^+1 for 
any value of j as 

7 ,-7,+i = 0(iV-§). (3.31) 

Thus we are able to prove that, for any polynomial type dispersion relation of the 
form (3.1), the asymptotic behaviours of 7, and 7, — 7 J+ i are given by Eqs. (3.29) and 
(3.31) respectively. Since these asymptotic forms do not depend at all on the parameters 
or coefficients present in the dispersion relation (3.1), they are indeed very universal in 
nature. Consequently, by applying the transfer matrix approach exactly like the cases 
of HS, PF and FI spin chains as have been discussed in Sec. 2, we can show that the level 
density for all ID vertex models and quantum spin chains associated with the polynomial 
type dispersion relation (3.1) asymptotically follow the Gaussian distribution (2.3) at 
iV — > 00 limit. 

It is worth noting that, as a by product of our analysis on the level density distribu- 
tion, we obtain Eq. (3.26) which gives a general expression for the order of a 2 for any 
dispersion relation of the form (3.1). In this context, it is interesting to ask whether 
such a formula exists for the asymptotic form of mean energy ji. Since the sign of the 
term £n{j) m Eq. (2.15) may change for different choice of j, such terms may cancel out 
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each other in the leading order. Due to this reason, it is difficult to get a general formula 
for the asymptotic form of p through the 'big 6' ordering of a function. Nevertheless, 
by using Eqs. (3.17) and (3.18), we derive a completely general formula for the upper 
limit of the asymptotic form of p as 



Let us now restrict our attention to those dispersion relations of the form (3.1), for 
which £n(J) ,s m Eq. (2.15) do not completely cancel out each other in the leading 
order. For such cases, the asymptotic form of p can be expressed through the more 
precise counterpart of Eq. (3.32) given by 



In particular, if the polynomial g x (p) in Eq. (3.15) does not vanish within the range 
< p < 1, then all £jvC?')'s have the same sign in the leading order and they can not 
cancel out each other. Therefore, the asymptotic form of p would be given by Eq. (3.33) 
for this case. By using Eqs. (1.4) and (3.15) we find that, for the cases of PF, HS and 
FI spin chains g x (p) is given by p, p(l — p) and p 2 respectively. Since these functions 
of p do not vanish within the range < p < 1, Eq. (3.33) is applicable for all of the 
corresponding spin chains. We have already noted that, x = 1 for the case of PF spin 
chain and x = 2 for the cases of HS and FI spin chains. Substituting these values of x to 
Eqs. (3.26) and (3.33), one can very easily reproduce the known results like a 2 = 6(iV 3 ), 
p = Q(N 2 ) for PF spin chain and a 2 = 6(iV 5 ), p = &(N 3 ) for HS and FI spin chains. 

4 Numerical results for a vertex model 

In the previous section we have analytically shown that the level density for a large class 
of ID vertex models and Yangian invariant quantum spin chains, associated with poly- 
nomial type dispersion relations of the form (3.1), asymptotically follow the Gaussian 
distribution (2.3). Here our aim is to compare this result with numerical calculations, by 
choosing some specific form of the dispersion relation (3.1) and taking different possible 
values for the related parameters. 

The numerical data for the level density distribution corresponding to any given 
dispersion relation can be obtained from the energy function E(s) defined in Eq. (1.6). 
In general, one can choose the spin path configuration s in many possible ways so that 
they lead to the same value of E(s), say, E{. Let us denote the number of such different 
spin path configurations by D^ rn \Ei). In the case of Yangian invariant quantum spin 
chains associated with the vertex models, this D (jn \E i ) represents the degeneracy factor 
of the energy level Ei. Choosing the values of N and m within some range, and using 



H = 0(N X+1 ). 



(3.32) 



n = e(N x+1 ). 



(3.33) 
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a symbolic software package like Mathematica, it is possible to obtain the values of 
all Ei and D^ m \Ei) for a given dispersion relation. With the help of such numerical 
data, one can draw the histogram for the corresponding level density and compare it 
with the Gaussian distribution. Furthermore, for the purpose of eliminating the effect 
of local fluctuations in the level density distribution, one may use the above mentioned 
numerical data to calculate the cumulative level density given by 



m 

E,<E 



and check whether this F(E) agrees well with the error function: 

E-/i 



C(E) = I e-^dx ' 



1 + erf 



(4.2) 



V2no* J-oo 2 
To get some specific dispersion relation with rather simple form, let us choose r = 1 
and s = 2 in Eq. (3.1), and set the values of three coefficients viz. / ,o, /i,o and /i,2 
equal to zero. Thus, we obtain 

£ N (j) = fo,d + fo,2j 2 + h,iNj , (4.3) 

where /o,i is an arbitrary real parameter, while /o,2 and /i,i are nonzero real parameters. 
Note that the degree of the polynomial £n{j) m Eq. (4.3) is given by x = 2. It may 
be observed that, by putting / 0;1 = 0, / ,2 = —fi,i = — 1 in Eq. (4.3), one obtains the 
dispersion relation for the HS spin chain. Thus the dispersion relation for the HS spin 
chain emerges as a special case of the dispersion relation (4.3). 

For the purpose of comparing level density with the Gaussian distribution (2.3), it is 
necessary to find out the values of corresponding mean energy and variance. Substituting 
£n{j) in (4.3) to Eqs. (2.15) and (2.16) respectively, we obtain explicit expressions for 
the mean energy and variance as 

A* = ^ (l " N(N - 1) [ (2/o,2 + 3A.0 N + (3/o,i - / 0)2 ) ] , (4.4a) 

(1 - -At )N(N- 1) (AN 3 + BN 2 + CN + D) , (4.4b) 
\ m 2 J 



a 2 = 



360 



where 



A = 6/ 2 2 + 5/ 1>1 (2/ 1)1 + 3/ , 2 ), B = 3/o,2(7/o,2 + 5/o,i) + 5/ lil (5/ 1 , 1 + 9/ , 2 + 4/ ,i), 
C = 5/0,1(2/0,1+9/0,2+ 10/i,i)- 49/ 2 )2 - 60/0,2/1,1, ,D = 5/0,1(5/0,!- 12/o,2)+ ll/oV 

Using Eq. (4.4b), we find that 
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It is interesting to note that the factor (6/q 2 + 10/i i + 15/o,2/i,i) appearing in the r.h.s. 
of the Eq. (4.5) can not vanish for any real nonzero values of the parameters /o,2 and 
fx t i. Thus from Eq. (4.5) it follows that cr 2 = 0(iV 5 ) for all allowed values of the related 
parameters, and this result is consistent with the general relation (3.26) for \ = 2- On 
the other hand, by using (4.4a), it is easy to see that the asymptotic form of \i depends 
heavily on the specific choice of the coefficients. For example, fi = 0(iV 3 ) for the case 
2/o,2 + 3/i 5 i 7^ 0, and the order of /i would be lower than N 3 otherwise. Again these 
observations are consistent with the general relation (3.32) for \ = 2. 

As a specific case, let us choose m = 2 and N = 40, and investigate the nature of level 
density distribution for the dispersion relation (4.3) by taking the related parameters as 
/o,i = 1, /o,2 = — 1, fi,i = 1- By using Mathematica, we get the numerical values of all 
Ei and D {m \Ei) for this case and plot a suitable histogram in Fig.l. From this figure 

0.001 



0.0008 
0.0006 

llT 

0.0004 
0.0002 


1000 2000 3000 4000 5000 6000 

Ei 

Figure 1: Continuous curve represents the Gaussian distribution G(E) and the his- 
togram represents the normalized level density distribution for the dispersion relation 
(4.3) with coefficients /o,i — 1, fo.2 = — 1, fi,i = 1, drawn for the case N = 40 and 
m — 2. 

it is evident that the histogram matches very well with the Gaussian distribution G(E) 
in Eq. (2.3), where the values of \i and a 1 are obtained by using Eqs. (4.4a) and (4.4b) 
respectively. Inserting the numerical values of D^(Ei) in Eq. (4.1), subsequently we 
calculate the cumulative level density distribution F(E) for the above mentioned case. 
In Fig.2 we plot F(E) at equidistant points and also the error function C(E) given 
in Eq. (4.2). Comparison from this figure shows that F(E) and C(E) are in excellent 
agreement with each other. We also calculate the mean square error (MSE) for this 
case and find it to be as low as 4.10 x 10~ 7 . In Table 1 we present the values of 
such MSE, which are calculated by taking different sets of values of / 0j i, /o,2, /i,i in the 
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1000 2000 3000 4000 5000 6000 

Figure 2: Continuous curve represents the error function C(E), while the circles rep- 
resent cumulative level density distribution F(E) for the dispersion relation (4.3) with 
coefficients /o,i = 1, /o,2 = — 1, /i,i = 1, drawn for the case TV = 40 and m = 2. 

dispersion relation (4.3) for a wide range of N (keeping the value of m unchanged). 
From this table we find that, for any particular choice of the three parameters, the MSE 



Sets of Parameters 


N = 20 




TV = 40 




TV = 50 




TV = 60 




TV = 70 




/o,i 


/o,2 


A,i 






















1 


1 


1 


8.17 x 10- 


-6 


1.18 x 10" 


-6 


6.56 x 10" 


-7 


4.08 x 10" 


-7 


2.74 x 10" 


7 


1 


-1 


1 


3.70 x 10" 


-6 


4.10 x 10" 


-7 


2.29 x 10" 


-7 


1.44 x 10" 


-7 


9.81 x 10" 


8 


1 


2 


1 


1.01 x 10" 


-5 


1.43 x 10- 


-6 


7.95 x 10" 


-7 


4.93 x 10" 


-7 


4.93 x 10" 


7 


1 


-1 


2 


2.46 x 10" 


-6 


4.09 x 10- 


-7 


2.33 x 10" 


-7 


1.48 x 10" 


-7 


1.00 x 10" 


7 


1 


-1 


3 


3.07 x 10" 


-6 


5.10 x 10- 


-7 


2.90 x 10" 


-7 


1.83 x 10" 


-7 


1.24 x 10" 


7 



Table 1: MSE for cumulative level density for the dispersion relation (4.3) with m = 2 

decreases steadily with the increase of TV. This fact clearly indicates that for m = 2 
the level density for the dispersion relation (4.3) asymptotically follows the Gaussian 
distribution. By using the same numerical method and taking different values of m, we 
have also studied the level densities for the dispersion relation (4.3) as well as a few 
other polynomial type dispersion relations and found that such level densities tend to 
the Gaussian distribution in all cases. 
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5 Concluding remarks 



By using the transfer matrix approach, we analytically show that the level density for 
a large class of ID vertex models and Y(sl m ) Yangian invariant quantum spin chains, 
associated with polynomial type dispersion relations of the form (3.1), asymptotically 
follow the Gaussian distribution. The universal nature of the asymptotic forms given in 
Eqs. (3.29) and (3.31), which do not depend on the parameters or coefficients present in 
the dispersion relation, play a key role in our analysis. We also present some numerical 
evidence in support of our analytical result. This result clearly implies that the Gaussian 
character of the level density distributions for these type of models is very robust in 
nature and probably originates from the underlying Yangian quantum group symmetry, 
rather than any specific choice of the dispersion relation. As a by product of our study 
on the level density distribution, we derive Eqs. (3.26) and (3.32) which give general 
expressions for the order of variance and mean of the energy function respectively, for 
any polynomial type dispersion relation. 

It is expected that the results obtained in this article would play an important role 
in investigating the density of spacings between consecutive energy levels of unfolded 
spectra for the class of Yangian invariant spin chains and ID vertex models associated 
with polynomial type dispersion relations. It should be noted that, for calculating the 
density of spacings, the energy levels are usually taken from the unfolded spectrum which 
has an approximately uniform level density distribution [46]. Such unfolded spectrum 
can be generated from the raw spectrum through a transformation involving the level 
density of the raw spectrum. Since we have shown that level density of all Yangian 
invariant spin chains and ID vertex models associated with polynomial type dispersion 
relations satisfy the Gaussian distribution, one can easily construct the corresponding 
unfolded spectra and analyze the density of level spacings for them [47]. It may also 
be interesting to study various thermodynamical properties and correlation functions 
for this type of Yangian invariant spin chains and vertex models through the transfer 
matrix approach. 
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